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SUMMARY 

Results of the stretch yo-yo feasibility and flight qualifica- 
tion tests are presented. These tests were conducted to prove 
the concept that the stretch yo-yo is a more accurate de-spin de- 
vice than the rigid yo-yo, and to verify the analytical development 
of the stretch yo-yo properties. Variations in the design param- 
eters and their effects on the final spin rate of the payload are 
noted in the analysis of the test results. The variables include 
initial spin rate, moment of inertia, and spring properties. A 
computer solution of the test payload equations of motion is in- 
cluded for comparison with the experimental results to confirm 
the mathematical analysis of the stretch yo-yo system. As a re- 
sult of the successful flight qualification tests a stretch yo-yo 
was flown on Ariel I (1962 o 1) in April 1962. 
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ANALYSIS OF THE DYNAMIC TESTS 
OF THE STRETCH YO-YO DE-SPIN SYSTEM 


by 

William R. Mentzer 

Goddard Space Flight Center 


INTRODUCTION 

The stretch yo-yo is a de-spin device that has the ability to compensate for errors in the initial 
spin rate and the moment of inertia of a payload. It is composed of an end mass, a helical spring, 
and, if necessary, a length of wire (Figures 1 and 2). The spring elongates under a load, giving a 
variable yo-yo length. The square of the yo-yo length varies directly as the initial spin rate, and the 
ratio of final spin to initial spin varies inversely as the square of the yo-yo length. This produces 
an essentially constant final spin rate. The concept of the stretch yo-yo was suggested by H. J. Cornille 
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Figure 1 —Stretch yo-yo consisting of wire, spring, and end mass. 
(The scale is in inches.) 
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Figure 2-'Stretch yo-yo consisting of a spring and end mass. (The scale is in inches.) 


(Reference 1). The analytical theory and design criterion were developed by Dr. J. V. Fedor (Refer- 
ences 2 and 3). 

Following dynamic analysis and development of design equations for the stretch yo-yo, the de- 
cision was made to fabricate and test stretch yo-yos on the Explorer XII (1961 vl) and Ariel I (1962 ol) 
type payloads. Three series of tests were conducted beginning with feasibility tests on the two pay- 
loads and concluding with flight qualification tests for the Ariel I payload. As a result of these tests, 
a stretch yo-yo was flown on the successful Ariel I spacecraft in April 1962. 

In the record of the experimental results, the effects of the several variables involved in stretch 
yo-yo design can be readily noted. These variables include optimum spring constant, preload, initial 
spin rate, spin moment of inertia, and material strength. 


OBJECTIVES OF THE TESTS 

The stretch yo-yo tests had several objectives which overlapped the three series of tests. The 
first series of tests was intended to demonstrate the feasibility of the stretch yo-yo de-spin system 
in compensating for errors in initial spin rate. This series also was used to examine the character- 
istics of a helical spring operating as a yo-yo. 

The second test series was planned as the flight qualification for the Ariel I stretch yo-yo made 
of National Standard Company's NS355 high-strength stainless steel. The yo-yos were to be tested 
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for sensitivity to changes in spin moment of inertia, and for performance at overspin conditions when 
the yo-yos would be operating near the yield point of the material. 

Because of an unexpected delay, the steel shipment did not arrive in time for the scheduled tests. 
Therefore, it was necessary to fabricate springs from conventional music wire. This change in spring 
material necessitated a change in the test series. The objects of the revised test series were to 
observe the effect of changes in moment of inertia, as originally planned, and to study the results of 
subjecting the spring to stresses greater than the yield point of the material. 

When the NS355 steel was received and the Ariel I springs had been fabricated, the third, and last, 
test series was conducted. These tests qualified the stretch yo-yo as Ariel I flight hardware. The 
yo-yos were tested at overspin conditions to see if they could withstand the loading. 


TEST APPARATUS 

Tests on the stretch yo-yo were performed in the vacuum facilities at Langley Research Center, 
Hampton, Virginia. The vacuum tanks were spheres 41 and 60 feet in diameter. The test structures 
were the portions of the outer shells of the payloads (Explorer XU and Ariel I) on which the yo-yos 
were positioned, and inertia plates for simulating the inertias of the complete payloads. A dc motor, 
sealed for vacuum operation, with an electromagnetic drive unit and an electromagnetic coupling and 
brake comprised the spin table drive system. The yo-yo firing signal was fed from an external 
manually operated circuit to the payload through a set of slip rings in the driven shaft of the spin 
table. Payload angular velocity was measured in the following manner. A disk with 32 equally spaced, 
radially protruding studs was mounted on the driven shaft; then, as the shaft turned, the studs 
generated pulses in an adjacent magnetic pickup. The output signal from this circuit was fed into a 
recording oscillograph. In order to obtain a complete time record of the operation, the declutching 
signal, firing signal, and yo-yo release signal also were fed into the oscillograph. 

The spin table and the mounting platform in the 60 foot vacuum chamber at Langley Research Cen- 
ter are shown in Figure 3. Figure 4 shows the Ariel I payload mounted on the spin table. The yo-yos 
can be seen in the picture. The details of the spin table with its drive mechanism and spin rate 
measuring device are shown in Figure 5. 

TEST PROCEDURE 

The original intent for the tests was to study only the ability of the stretch yo-yo to compensate 
for errors in the initial spin rate of the payload and errors in payload moment of inertia. But, as a 
result of yo-yo fabrication problems, the experimenters were forced to consider the effects on the 
final spin rate of preload in the yo-yos, of deviation from the optimum spring constant, and of spring 
loading in the region of the elastic limit of the material. Because of the increased complexity of the 
program, it was necessary to be careful in pairing the springs for the test. The spring scale and the 
preload were determined for each spring to be tested. Heat treating and flexing were performed on the 
springs in an attempt to decrease the preloads. One spring from each group fabricated was statically 
loaded until permanent deformation occurred. This was done in order to obtain a strength limit to 
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Figure 3— The Ariel I payload mounted on the spin 
table in the 60 foot vacuum chamber. 


compare with the expected maximum yo-yo 
tension in the upcoming tests. With regard 
to the two factors affecting the tests, spring 
constant and preload, a deviation from the 
optimum spring constant was the variation 
most critical to the test results. In view of 
this fact, the springs with essentially equal 
spring constants were paired for testing. 
The pairs were then tested in the order of 
smallest to largest preload. 

After the yo-yo springs were paired, the 
average values of spring scale and preload 
for a set were used to compute the end mass 
(Reference 2). The yo-yos were assembled 
so that both springs in each test were of the 
same length. 

The tests were conducted at an absolute 
pressure of 10 mm Hg to minimize the at- 
mospheric drag effects on the yo-yos. Pay- 
load spin rate was determined from the re- 
cording oscillograph records (Figure 6). 



Figure 4— Closeup of the Ariel I payload showing 
the yo-yos in place. Weights at the base of the 
structure are for simulating the entire payload 
moment of inertia. 


MAGNETIC CONTACT 
FOR SPIN MEASUREMENT 



Figure 5— Details of the spin table. 
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FIRING SIGNAL 



The oscillograph output consisted of a plot of pulses from the magnetic contact on a known time scale. 
Also included on the record were the declutching signal, the yo-yo firing signal, and the yo-yo release 
signal. 

The spin table was constructed so that one pulse from the magnetic contact took 1/32 of a revolu- 
tion of the payload and that the time scale of the recorder chart r was 1/100 of a second per division. 
The spin rate at the midpoint of the time interval between any two pulses was found from the relation 

N 

W = 1.875 T , 


where 


w = spin rate in rpm, 

N = number of pulses, 
t = total time between pulses in seconds. 

DISCUSSION OF TEST RESULTS 

The stretch yo-yo tests were considered successful from several standpoints. The concept of 
the stretch yo-yo as a more accurate de-spin device than the rigid yo-yo was proven, and the analyti- 
cally developed design criteria were verified (Reference 1). The difficulties of operating with a non- 
optimum yo-yo system were evaluated. Finally, the yo-yo was qualified as flight hardware for Ariel I. 

The results of the stretch yo-yo tests are recorded in Table 1. The tests are grouped by series 
and by payload and are numbered consecutively from the first successful test on each payload. Sys- 
tem design parameters and actual test conditions have been tabulated in order that the effect on final 
spin rate of variations of test conditions from design values can be readily noted. 

The Explorer XII type optimum spring scale tests clearly demonstrated the ability of the stretch 
yo-yo to compensate for errors in initial spin rate. In this series of tests, spin- up errors of ±20 
percent were reduced to within ±1.5 percent of the design final spin rate. The fine performance was 
attributed to the fact that the yo-yo springs had the optimum spring constant and no preload. The 
yo-yos used in these tests were of the type shown in Figure 1. 

The importance of the optimum spring constant in the proper functioning of the stretch yo-yo 
de-spin system was illustrated by the last three tests on the Explorer XII configuration. These tests 
were necessitated by the manufacturing of preloaded springs with incorrect spring constants. The 
tests were conducted at a higher spin level than tests II through IV because of apprehension concern- 
ing the functioning of a preloaded yo-yo. Spin limits were selected and calculations were made by using 
the measured values of preload and spring constant, which were nonoptimum for the spin rates used. 

The percent error shown for these tests is large when referenced to the 30 rpm design final spin. 
If the results are referenced to test V, with a 26.8 rpm final spin rate, as design conditions, then the 


Table 1 

Stretch Yo-Yo Test Summary. 
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*One yo-yo had zero preload. 
tYielding occurred. 




overspin error is + 2.7 percent. The underspin error, -23.6 percent, is still large but is attributed 
to the fact that the test was made with yo-yos of unequal preloads, differing by 2.6 lb. 

Ariel I yo-yo tests II and III proved the feasibility of the one-half wrap stretch yo-yo system 
(Figure 1) and led to the use of the stretch yo-yo on this satellite. Tests IV and VII demonstrated the 
repeatability of the yo-yo results and, along with test VI, gave a picture of the system operating at 
the design moment of inertia. 

The design final spin rate for these tests was computed from equations that included the effects 
of preload in the yo-yo springs. The value of the preload used in calculations was that of the springs 
used in test IV. Test results agreed closely with theoretically predicted final spin rates, the deviation 
resulting from the fact that the spring constants varied slightly from the optimum spring constant. 

In tests V and XI the springs yielded because the load exceeded the yield point of the material. 
The yielding produced a spin rate that was much lower than theoretically predicted. This situation 
was permitted to occur since the springs being tested were of lower strength limits than flight hard- 
ware. The experimenters were afforded the opportunity to observe the behavior of the de-spin system 
when yielding occurred and the design equations no longer applied. 

Tests VIII through XII verified the theoretical prediction that the stretch yo-yo would be relatively 
insensitive to variations in spin moment of inertia when compared with the behavior of the rigid 
yo-yo. Moment- of- inertia values of ± 20 percent of the design spin axis moment of inertia were used 
in the tests. 

The stretch yo-yo tests were concluded with the Ariel I flight qualification tests, XIII through XV. 
In these tests, the NS355 high-strength stainless steel springs were used. The tests were conducted 
at design initial spin rate and then at + 15 and + 20 percent of design initial spin rate to determine 
whether the yo-yo could withstand the high loading at overspin conditions. Test XIII, with nominal 
initial spin rate, gave a final spin slightly higher than design final spin. In the overspin tests, the 
final spin rates were very close to design values. As a result of this test series, the stretch yo-yo of 
high-strength stainless steel was accepted as flight hardware for Ariel I which was launched suc- 
cessfully in April 1962. 

An analytical method for determining the actual final spin rate of a satellite is given in Refer- 
ence 2. The data necessary for performing the calculations include spring properties, design values 
of spin rate and inertia, and actual values of spin rate and inertia. 

The equations of motion of the test payload were programmed for solution on an IBM 7090 digital 
computer as a verification of the test results. Correspondence between this solution and the test re- 
sults proved the validity of these equations of motion. The equations were developed for phase 1 of 
the yo-yo operation, when the yo-yo is unwinding and is tangential to the payload. Spin table friction 
was included in the analysis. The development of the equations and the computer program is included 
as Appendix A. 

A graph of the computer solution of phase 1 for one test is included with the plot of the entire de- 
spin test in Figure 7. The plot shows the angular displacement of the payload, in radians, versus time. 
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The test data points are the pulses from 
the magnetic contact that were recorded 
by the recording oscillograph. Thirty-two 
successive points define one revolution 
of the payload. The end of phase 1 is de- 
termined by the computer program, and 
the end of phase 2 is indicated on data 
tape by the yo-yo release signal. At the 
end of phase 1 the computed angular dis- 
placement for test XIII of February 7, 
1962, exceeded the recorded value by 
2.42 percent. 

Phase 2, when the yo-yo moves from 
a tangential to a radial position, was not 
analyzed because of algebraic complexity 
and similarity to phase 1 analysis. It was 
felt that the close agreement between the 
computer solution of the phase 1 analysis 
and the test results, combined with the 
similarity of the methods of solving phase 
1 and 2 equations on a computer, justified 
the omission of the phase 2 analysis from 
this report. 



Figure 7— Angular displacement vs. time for Ariel I 
test XIII, February 7, 1962. 
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Appendix A 


Computer Solution of Phase 1 Equations of Motion 


The theoretical development and the design criteria for the stretch yo-yo de-spin mechanism 
functioning on an orbiting spacecraft have been developed by J. V. Fedor.* The stretch yo-yo tests 
were performed to demonstrate the feasibility of the stretch yo-yo system and to verify the design 
parameters that had been developed. Since the observed final spin rates of the payloads corresponded 
to the theoretical predictions, the tests were considered successful. 

The differential equations of motion of the test payload on the spin table were solved on an IBM 
7090 computer for comparison with the experimental results to verify the analysis of the stretch yo-yo 
system. The equations were developed for phase 1 of the yo-yo operation when the yo-yo is unwinding 
and tangential to the payload (Figure Al). In the analysis, the effects of gravity are considered negli- 
gible because of the short operating time of the yo-yo, 0.50 to 0.75 seconds, and because of the exist- 
ence of a component of the tension in the yo-yo spring which opposes the gravitational force. The 
validity of this assumption has been demonstrated 
in previous rigid yo-yo de-spin tests where test 
results agreed closely with theoretical calcula- 
tions in which gravitational effects were neg- 
lected. Tests were conducted at a pressure of 
10 mm Hg, which corresponded to an altitude of 
100,000 feet and made the effect of atmospheric 
drag negligible. 

The coordinate system for phase 1 analysis 
is shown in Figure Al. Because of the symmetry 
of the system only one yo-yo is shown. 

The stretch yo-yo de-spin system behaves 
as a rigid yo-yo until the preload is overcome by 
the force in the yo-yo. During this phase of op- 
eration the Lagrangian L , which is the kinetic 
energy of the system, is 1 

L = y 14> 2 + ym (l 2 e 2 + a 2 <*> 2 ) , (Al) 


y 



Figure Al— Phase I coordinate system. 


*Fedor, J. V., "Analytical Theory of the Stretch Yo-Yo for De-Spin of Satellites/’ NASA Technical Note D-1676, April 1963. 

^ Fedor, J. V., "Theory and Design Curves for a Yo-Yo De-Spin Mechanism for Satellites/ 1 NASA Technical Note D-708, August 1961. 
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where 


I = spin moment of inertia less yo-yo masses, 

<t> - spin rate of payload, 

6 - generalized coordinate, 
m = mass of both yo-yos, 
a = radius of yo-yo fixture, 
l = length of yo-yo unwound. 

The length of wire unwound at any time in rigid yo-yo operation is 


l = a(0 ~ 0) . (A2) 

The mass of the yo-yo system, including the spring mass, is given by an approximate equation that 
accounts for the distributed spring weight, 

m = m o + nr - (A3) 


in which (pa/3) (0 - <p) is the mass of yo-yo spring unwound for both yo-yos. The quantity p is twice 
the mass density of one helical spring and m 0 is the total end mass of the yo-yos. The equations of 
motion in terms of the Lagrangian are 

Q 4> ’ (A4) 

Qe • (A5) 

Evaluation of these expressions, in which the generalized forces and Q d consist only of the fric- 
tion torque in the spin table, -si andQ 5 = 0, yields 

■Jj: [^(i + ma 2 ) + ma IB 2 - ~sl , (A6) 

3f(m l 2 6) - ma IQ 2 = 0 , (A7) 


d. 

> <JL\ 

dl 

■atl 


~ d<f> 

d i 


dl 

dtl 

\db) 

~ d0 


where s is the deceleration rate of the payload due to friction measured at t =0. 

In order to obtain a difference equation for ^asa function of time, integrate Equation A 6: 

r k+i 


£(l +ma 2 ) <)>J - - f (si +ma Id 1 ) dt . 


Then 


l 1 +m k + 1 a 2 ) 0 k + 1 - (l+m k a 2 ) 4 - sIAt- m k aZ k 0 k 2 At . ( A8 ) 

Represent 4> k +1 by a forward difference, (<?b k+2 ~^ k+ 1 )/At, and then solve Equation A8 for<£ k+2 : 

At 


^k +2 “ ^k + l 


1 +m k+l 32 


[(! +m k a 2 ) 4 " si At - m k a *k 4 2At 


] 


(A9) 
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Performing the indicated differentiation of the second equation of motion yields 

m'S + mb = T (ma 9 2 ~ 2m Id) , (A10) 

and the left-hand side of this expression is an exact differential which can be expressed as 

^ lm£» = x (ad - 21) . (All) 

Integration of Equation All, by the method of Equation A8, combined with the replacement of 6 k + l by 
a forward difference, gives the relation for £ k+2 : 


c + 2 ^k + i + ^k m k + 1 


1 + (2a <j> k - a£ k ) 


At 


kJ l 


kJ 


(A12) 


The computer solution is started by calculating the values of $ and 6 for t = 0. Att -o,i £ o 
and the right-hand side of Equation A12 becomes infinite; thus, an expression must be developed to 
evaluate 0 k+2 at the time t = o. 

From the expanded form of the second equation of motion, Equation A10, an expression for e 0 is 
determined. Substitution of initial conditions into Equation A10 yields 


= 24>o • 


(A13) 


Application of L'Hospital's rule to the term (2a <£ k - a9 k ) At/l k from Equation A12 yields, at t =0, 

<P 0 

From the differentiation in Equation A6 the general expression for the payload angular acceleration is 


-m k a 2 - m k a £ k 2 - si 
I + m k a 2 


(A14) 


which can be evaluated at t - 0: 


- f tj - si 

I + m Q a 2 


(A15) 
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Upon differentiating in Equation A7 and solving for 0, it is found that this expression is undefined at 
t = 0; and L'HospitaUs rule must be applied again. The resulting expression for 0 at t =0 is 


e 


0 


‘ 2 P ai 0 2 4<ji> 0 

9% + T ' 


(A16) 


A combination of ( 24> 0 - 6 0 ) At /4> 0 with Equation A12 yields an expression for d 2 in terms of the ini- 
tial conditions: 


where 6 i = & g At + & 0 . 


0i + 


l +- 


At(2^ o -0 o )" 


nij e o ^ 


(A17) 


The total force in the yo-yo system during the time when it performs essentially as a rigid yo-yo 
is given by 


F k - m k( a ^k + K K 2 ) + K a A • ( A18 ) 

The rigid yo-yo operation is terminated when F becomes greater than 2F 0 , the total preload in the 
yo-yo springs. 

After the tensile force in the yo-yo has exceeded the preload, the yo-yo functions as a stretch 
yo-yo. The Lagrangian for the stretch yo-yo in phase 1 of yo-yo operation is* 

L = \ 102 + i m(z 2 e 2 + a 2 <£ 2 + S 2 - 2a0s) - (kS 2 + 2F Q s) , (A 19 ) 


in which 

B = deflection of a spring, 
k = spring constant, 

F 0 = preload in one spring. 

The length of the stretch yo-yo at any time t is a function of the deflection and the amount of spring 
unwound; thus 


l = a(B-4>) + S ; (A20) 

whereas the mass remains dependent only on the length of wire unwound, 

pa 

m - m 0 + “3“ {6 -<t>) . 


"Fedor, J. V., "Analytical Theory of the Stretch Yo-Yo for De-Spin of Satellites," NASA Technical Note D-1676, April 1963. 
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Writing the equations of motion explicitly from the Lagrangian results in: 


^ [(i +ma 2 )^-ma §] +m a W 2 - - si , (A21) 

l 2 0) - maid 2 = 0 , (A22) 

^(m S - m 16 2 + 2kS + 2F Q = 0 . (A23) 

To determine <j> as a function of time for stretch yo-yo operation, integrate Equations A21 and A23 in 
the manner of the preceding work to get 

(i + m k+1 a 2 ) 4> k + 1 - m k + 1 a§ k + 1 - (i + m k a 2 ) <£ k + m k a§ k = “ m k a l k 6 2 At - si At , (A24) 

m k + i ^k + i “ m k + i a< ^k + i “ m k ^k + m k ac ^k " m k ^k ^k 2At “ 2kS k At - 2F 0 At (A25) 

Solve Equation A25 for m k + 1 S k + 1 and substitute into Equation A24, simultaneously replacing <£ k + 1 by a 
forward difference, (<£ k + 2 -<£ k + 1 )/At , to get an expression for the angular coordinate of the payload: 

A + 2 = ^k + 1 + X ( l A - si At - 2ak8 k At - 2aF 0 At) . (A26) 

A treatment of Equation A22 by the same method of integration and substitution yields 

= g k+. + m At , 2 Ka; k V At+ % V K) . /A ,,v 

where £> k + 1 = (e k + 2 - 0 k + l )/At as above. 

To determine the difference equation for S, differentiate in Equation A23 as indicated and then 
simplify to the form: 

d • • 

(mS) - ma<j> + ma 4> + m W 2 - 2kS - 2F 0 (A28) 

Integrate this expression and substitute a backward difference for 8 to get 

S k + i = S k + li^[ m k S k +(m k 4 k +m k a^ k +m k Z k 6> k 2 - 2kS k -2F 0 ) At] . ( A2 g) 
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This relationship is dependent upon which is determined by differentiating in Equation A21 as indi- 
cated and substituting from Equation A28: 


= - y( sI + 2ka5 k + 2F o a ) • (A3Q) 

If there is no preload in the yo-yo springs, the device will function as a stretch yo-yo from the start 
of the de-spin operation. For this reason it is necessary to be able to evaluate the expressions for 4>> 
0 , 8, and '<p at t = 0 . Differentiation in Equation A22 gives 

mO + m& - ~j{mad 2 ~ 2mi&) , (A31) 


and substitution of initial conditions yields 


' e o = 2 ^o ■ 

which is also true for rigid yo-yo operation. At t = 0, Equation A30 reduces to 


A) " " s > (A3 2) 

which is compatible with the physics of the problem in that the payload is decelerating because of 
friction at t = o. The deflection equation becomes 

At 2 / . . ■■ . 

b i ~ mj K at ^o +m o a< t>o) ’ (A33) 

at t = 0 . Equation A26 for the payload angular coordinate takes the form: 

^2 = <*>1 + T^(l^o " slAt ) . (A34) 

for t = 0. Equation A27 for 9 cannot be evaluated at t = 0; thus, another formulation of the problem 
must be considered. Rewrite Equation A31: 


(m<9) - y(ma^ 2 ~ 2mld) 

Then integrate and make a forward difference substitution for e: 


k + 2 


c + 1 + ^k At 


1 + (2a^ k -a0 k -2S k )f^ 


(A3 5) 
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This equation cannot be directly evaluated at t = 0 . Apply L'Hospital's rule and initial conditions to 
the quantity (At/Z k ) (2a <£ k - ad k - 28 k ) to get (2<ji 0 - £ 0 ) At /<p 0 , where 9 0 is found from Equation A31 by 
the same method used in deriving Equation A16: 


- 2 pa <£ 0 2 4 4> 0 

+ ~ 


And we have for Equation A3 5 


+ 




¥7^0 At 


(A36) 


(A37) 


with (p 0 from Equation A32 and 0 0 from Equation A36. 

Phase lof the yo-yo operation is terminated when the entire yo-yo is unwound from the payload. 
The computer terminates the program when the relation 


[ B “ a K -0k)] - 0 (A3 8) 

is satisfied, where B is the length of the yo-yo and a [o k -<£ k ) is the length of yo-yo unwound at time 
t = V 

An analysis of phase 2 of the yo-yo has been omitted from this report. Since the equations are 
similar to those for phase 1* and the computer solution of phase 1 was sufficiently close to the test 
data, it was felt that solution of the phase 2 equations would not add sufficient information to the re- 
port to warrant the additional effort involved. 

The computer program is outlined on the following pages in the form of a block diagram, Figure 
A2, and a listing of the program in Fortran for the IBM 7090, Table Al. The nomenclature for the 
program is as follows: 

<t> = PHI, 

6 = THETA, 

S = DELTA, 
l = EL, 
m = EM, 

F = FORCE, 
t = TIME, 


and in the input data: 



a = A, 

1 — R 

^ SPRING 9 

(i + mg a 2 ) 4> 0 - PI3, 

0 O = OMEGAO, 

At = DTIME, 

si = PI4, 

m 0 = EMO, 

i = pii, 

2k = PI5, 

2F 0 = P0, 

pa/3 = PI2, 

2F 0 = PI6. 


•Fedor, J. V., "Analytical Theory of the Stretch Yo-Yo for De-Spin of Satellites,” NASA Technical Note D-1676, April 1963. 


OMEGA, 

<£ = ALPHA. 

ETA, 

6 = BETA, 

ZETA, 


DEEM, 
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READ IN 



A, OMEGAO, EM0 


WRITE OUT 

P0,B,DTIME 
PI1— PI6 


(same ) 


j£ 


COMPUTE 
TIME = 0.0 
PH 12, OMEGA1 
THETA2, ETA1 
XI, ELI, EMI 



COMPUTE 

TIME 

EM2, EL2, X2 
PHI3 

OMEGA2 

SMALR2 



COMPUTE 

ALPHAO 

BETAO 

THETA3 



PHI2 = PHI3 


OMEGA! - OMEGA2 
ALPHA! = ALPHA2 


THETA2 = THETA3 


ETA! - ETA2 


EL! = EL2 


EM! = EM2 


XI =X2 




COMPUTE 

ETA2 

ALPHA2 

FORCE2 


1 

l 

’20 

WRITE 

(calcul 

1 

OUT 

ations ) 

1 


WRITE OUT 

(end of rigid yo-yo ) 


PHI! = PHI2 


PHI2 = PHI3 


OMEGA! = OMEGA2 


ALPHA! = ALPHA2 


THETA! = THETA2 


THETA2 = THETA3 


ETA1 = ETA2 


EL! = EL2 


EMI = EM2 


DELTA! = 0.0 


ZETA1 =0.0 


no 


41 


COMPUTE 
TIME 
EM2, EL2 
PH 1 3, OMEGA2 
SMALR2, ALPHAO 
DEEM! 
DELTA2 
ZETA1 
THETA3 


COMPUTE 
TIME, EM2 
DEEM1 

PHI3, OMEGA2 
SMALR2, DELTA2 
ZETA2, EL2 
ALPHA2, THETA3 



Figure A2— Flow chart for the solution of the phase 1 equations of motion of the stretch 
yo-yo for the IBM 7090 digital computer. 
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Table A1 


Fortran Computer Program for Phase 1 . 


* W.R.MENTZER MECHANICAL SYSTEMS BRANCH G.S.F.C. 

* XEO 

-* C A R PS COLUM N 

* LIST 8 

C ST R ETCH TO YO PHASE 1 

1001 FORMAT ( 3XF6.4*2XF8.4»3XF8.6*3XF7.4*3XF7.4*3XF7.5*2XF6. 3/4XF9.8*3XF 

C 9 .6*3XF10.8«3XF 7 .3»3XF6.31 

1002 FORMAT(6X1HA.5X6HOMEGA0.6X3HEM0.9X2HP0.7X1HB*8X5HDTIME*6X3HPI1*6X3 
C H P 1 2 i9X3HP I 3»9X3HP 1 4 »6X ? H P I 5 *6X3HP Ife/3XEfe«4t2XF8«4»3XFB.6* ' 3XF 7 .A ' »3 

CXF7.4.3XF7.5*2XF6.3«2XF9.8*3XF9.6.3XF10.8*3XF5.2*2XF&.3///) 

1003 FORMAT < 1H1 *20X5 1H P ROGRAM F OR T H E TE ST O F T H E S T RE I' CH YO VO PHASE 

Cl //5X4HT IME*4X10HSPIN RAT 1 0 , 4X9HSP I N RATE *6X 3HPH I , 6X 1 2HACCELERA T I 

CON*5X5HTHETA*8X6HLENGTH*5X10HDEFLECTION*5X4HMA55»6X5HFORCE/5X5HT5F 

CC> *17X9H< RAD/SEC) . 5X5H( RAD >» 5X1 3H { RAD/SEC/SEC ) *4X5H( RAD) « 9X4H C FT j * 
C9X4H t F T ) * 7 X 7 HT ' S C UGS ) *5X4H< LB)/// ] 

1004 FORMAT ( 3XF7,4*2XF11«8*3XF10.6*3XF10.7»3XF10.7*3Xf10.7*3XF10.6*2XF1 

Clt 8.3X r9 »8»2X F7 .4/) — 

1005 FORMAT ( ///6X31HEND OF OPERATION AS RIGID YO YO ) 

1006 FORMAT! 3XF7.4.2XFirrr.3XF10.6.3XF10.7.3XF10.7.3XFlO'r7*3'XFTo;'6"*2X'FT 
C1.8»3XF9.8/) * 

1007 FORMAT (///6XT4HEND OF PHASE 1 1 

1008 FORMAT ( 6X12HNEGAT IVE ELI) 

RSAD TNPttT rAPE 2 » 1001 .A.OMEGAO , EMO *P0*B »DT IME . PI 1 »P I 2 »PI 3 *PI4 t P 

CI5.PI6 

WRITE OUTPUT TAPE 3* 1002 ♦ A*OMEGAO * EMO »P0 *B *DT IME »P1T *P I 2 »PI 3 »P 14 *P 

CI5.PI6 

y inTE O Cf T PU T T APE 3 >1003 

J = 1 

1 - t IME= 0.0 

FORCE1=0.0 

EMl-stMO 

EL1=0.0 

O M EGA 1*O M EO AO 

PHI 1=0*0 

— PM I2=0MEGA1»DT IME 

ETA1 = 2 • 0*OMEGA0 

TMETAl^O.O - 

THETA2=E T A1»DT IME 

IF (PO) 9*110*10 

9~CALL'FXn 

110 TIME=TIME+DTIME 

DELTA1 = 0.0 

ZETA1=0.0 

EM2 a EM0+P 1 2* ( THET A2-PHT 2 ) 

EL2=A*< THETA2-PHI2 ) 

PHI 3*PHI Z+ ( DT IME/FT1 1 STFTT*0MEGA1-P l 4*DTTME ) 

OMEGA2 = (PHI3-PHI2 ) /DTI ME 

SMALR2=OMEGA2/OMEGA'0 “ 

ALPHA1= -P 14/P 1 1 

DEE Mls C E M2- E M1 ) /P TT M E 

DELTA2* (DTIME»*2/EM2 )*(DEEM1*A*0MEGA1+EM1*A*ALPHA1 ) 

Z ETA2= ( DELTA2-DELTAn/DTl ME ~ ~ — 

THETA3=THETA2+DTIME*ETA1*(EM1/EM2)»(1.0+(DTIME/OMEGA1 )*(2.0*ALPHA1 

— — C-BETA1 ) I 

ALPHA2=-( l.O/PI 1)*(PI4+A*PI5*DELTA2+PI6*A) 

GO 'TO" 50 
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Table A1 (Continued) 


10 TIME=TIME+DTIME 

EH2.=EM0.*PT2*< THETA 2-PHI 24 

EL2=A*( THETA2-PHI2 ) 

PHI3=PHI2+(DTIME/X2)*<X1*0MEGA1-PI4*DTIME-EM1*A*EL1*ETA1**2*DTIME ) 

Q MEG A 2 = ( PH l-^-P-Ht 2T/-DT 1 ME - 

SMALR2=OMEGA2/OMEGAO 

TF - t- ELl ) 13.14.15 - - 

13 WRITE OUTPUT TAPE 3*1008 

CAU.--EXTJ. - 

14 ALPHA1=-(PI 2*A**2#OMEGA 1**2 +PI41/X1 

' BETA Q#PI-?*C1MEGA 1##21 / <-3.G*EMO \* (4.0* ALPHA 1 1/3.0 

THETA3=THETA2+DTIME*ETA1*( EM1/EM2 ) *( 1 .0+ ( DT IME/0MEGA1 ) * ( 2 . 0*ALPHA 1 

C -B.ETA1 1 1 - — - 

GO TO 16 

15 lHElAl=-I.hi£JA23-0IlM£*E TA1#1EM1/EM2 1* (1.0+12. 0* A* OMEGA 1-A*FTA1 ) * ( DT I 

CME/EL1 ) ) 

GO TO 16 . 

16 ETA2=(THETA3-THETA2)/DTIME 

ALPHA2=. .-.( P.I2*A**2*OMEGA2* ( ETA2-OMEGA2 ) +EM2*A*EL2*ET A2**2+P 1 4 ) /X2 

EM3=EM0+PI2*( THETA3-PHI 3) 

...Q££M2= .( EM3 -£M2 1 /0.T.I ME ... . . 

FORCE2 = EM2*< A*ALPHA2 + EL2*ETA2**2 ) + A*DEEM2*OMEG A2 

2.0. IF.L5L - J) 21.22*23 . 

21 CALL EXIT 

.22 ..WRITE OUTPUT TAPE 3.1003 

WRITE OUTPUT TAPE 3*1004. T I ME .SMALR2 .0MEGA2 * PH 1 2 » ALPHA2 ♦ THE T A2 » FL 

C2..QE1IA2 *£M2 j FORCE2 . 

J = 1 

GO TO 25 - 

23 WRITE OUTPUT TAPE 3*1004. T I ME .SMALR2 .0MEGA2 . PH 1 2 * ALPHA2 . THE TA2 * EL 

C-2-*£>E-LTA2»EM2 .FORCE 2 

J = J + 1 

60 - TO— 2 5— - - 

25 I F ( PO-FORCE2 ) 40.40.30 

30- P++l-2=PH I3 - ..... 

0MEGA1 =0MEGA2 

ALPHA 1=ALPHA2 ....... 

THET A2 = THETA3 

^1AT=€-TA2- 

EL1=EL2 

EKa=EM2 

X1 = X2 

GO TO 10 . 

40 WRITE OUTPUT TAPE 3.1005 

PH I 1 =PH 1 2 

PHI 2 = PH I 3 

OMEGA 1 = 0MEGA2 - - .. 

ALPHA1=ALPHA2 

THE T A1 s.IHE T A2- - — - 

THETA2=THETA3 

£TA1=EJ.A2 - - - 

EL1=EL2 

EHls£M2 - - - 

DELTA1=0.0 

2EIAl-=0-0 - 

WRITE OUTPUT TAPE 3.1003 

41 TIME=TIME+OTIME 

- EM2— EMO-fcR J-2-*tTHFTA2-^PH I 2 )- - - - 

DEEMl=(EM2-EMi l/DTIME 
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Table A1 (Continued) 


PH I 3=PH 1 2 + (DTI ME/P II ) * ( P I l*OMEGAl-P I4*DT I ME~A*P 1 5*DELTA1*DT I ME-P 1 6 

- C*A»QT IM E) . - _ 

0MEGA2= ( PHI 3-PHI 2 ) /DTI ME 

SMAI R2sDMFfiA2/OMFGAO 

DELTA2=DELTai+( DTIME/EM2 )*<EM1*ZETA1+(DEEM1*A*0MEGA1+EM1*A*ALPHA1+ 

C E M 1 *E 1 1 »E T A- 1 * * g -~p.j-5A-DE-4.-T A 1 -P 1-6- )»0 TIME ) - 

ZETA2= (DELTA2-DELTA1 ) /DTI ME 

EL2 —A * ( - THE 1A 2 - PH1 2 1+ DELTA-2 — 

ALPHA2= -( 1.0/PI 1 )*<PI4+A*PI5*DELTA2+PI6*A) 

I H ETA3-THET - A? * ( -l>T4-ME-/-( &M? »EL-2*»?-) I'M EMT***EbT*ETAl**2*DT IME+EM1-»EL 

C1**2*ETA1 ) 

E - T -Ar2= ( T H E -T A 3 - T H E 7-A2 F/OT-LME- - - 

5C I F ( 51 - J) 51*52*53 

54- CALL EXIT 

52 WRITE OUTPUT TAPE 3*1003 

C2»DELTA2*EM2 


GO TO 60 


-4*-1006*- 4-IME-*SMALR2 *0MEGA2 *PHI 2* ALPHA2 * THETA2 *EL 


C2 *DElTA2 *EM2 


l=J+J 

GO TO 60 

6Q— LF.XB-tA* ( THETA2«-PH 1244 -90*90.40 — 

80 PHI 2 = PH I 3 

OMEGA 1=0MEGA2 

ALPHA1=ALPHA2 

TH£IA2.?..IHE1A3 

ETA1=ETA2 

FI 1=EL2 

EM1=EM2 

DELTA1 = DELI-A2 

ZETA1=ZETA2 

GO— T O AX 

90 WRITE OUTPUT TAPE 3*1007 

CAL L E X IT 

END 

Jt— DATA 




THE DATA CARDS FOR THIS PROGRAM USING THE DATA FROM TEST III OF 
FFR. 7 1962 ABF...A£ . FOLLOWS 

A -O MEG A 0 EMO PO B DT IMF PIT — 

.942 16.8300 .033850 3.8000 2.3650 .0002 2.885 


P 1 2 P 1 3 P 14 P 1 5 

Milan ASL..Q60073 *-40181584 -2 3*460- 


P 1 6 

3*800 


NASA- Langley, 1963 


G-394 



